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It is well known that the KdV equation has an infinite set of conserved quantities. The first three 
are often considered to represent mass, momentum and energy. Here we try to answer the question 
of how this comes about, and also how these KdV quantities relate to those of the Euler shallow 
water equations. Here Luke’s Lagrangian is helpful. We also consider higher order extensions of 
KdV. Though in general not integrable, in some sense they are almost so, these with the accuracy 
of the expansion. 
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I. INTRODUCTION 

There exists a vast number of papers dealing with the 
shallow water problem. Aspects of the propagation of 
weakly nonlinear, dispersive waves are still beeing stud¬ 
ied. Last year we published two articles a a in which 
Korteveg-de Vries type equations were derived in weakly 
nonlinear, dispersive and long wavelength limit. The sec¬ 
ond order KdV type equation was derived. The second 
order KdV equation [3, 3 ], sometimes called "extended 
KdV equation", was obtained for the case with a flat 
bottom. In derivation of the new equation we adapted 
the method described in [4]- In [2j, an analytic solution 
of this equation in the form of a particular soliton was 
found, as well. 

It is well known, see, e.g. @-(1], that for the KdV equa¬ 
tion there exists an infinite number of invariants, that is, 
integrals over space of functions of the wave profile and 
its derivatives, which are constants in time. Looking for 
analogous invariants for the second order KdV equation 
we met with some problems even for the standard KdV 
equation (which is first order in small parameters). This 
problem appears when energy conservation is considered. 

In this paper we reconsider invariants of the KdV equa¬ 
tion and formulas for the total energy in several different 
approaches and different frames of reference (fixed and 
moving ones). We find that the invariant I^\ some¬ 
times called the energy invariant, does not always have 
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that interpretation. We also give a proof that for the sec¬ 
ond order KdV equation, obtained in [lHl, ./A) r i 2 dx is 
not an invariant of motion. 

There are many papers considering higher-order KdV 
type equations. Amongthem we would like to point out 
works of Byatt-Smith Q, Kichenassamy and Olver Ea, 
Marchant a [rnii, Zou and Su EH; Tzirtzilakis et.al. 
0 and Burde E3- It was shown that if some coef¬ 
ficients of the second order equation for shallow water 
problem 0 are diferent or zero then there exists a hier¬ 
archy of solition solutions. Kichenassamy and Olver 0 
even claimed that for second order KdV equation solitary 
solutions of appropriate form can not exist. This claim 
was falsified in our paper Q where the analytic solution 
of the second order KdV equation 0 was found. Con¬ 
cerning the energy conservation there are indications that 
collisions of solitons EM! which are solutions of hig her 
order equations of KdV type can be inelastic [HI Ea |. 

The paper is organized as follows. In Section ITIl several 
frequently used forms of KdV equations are recalled with 
particular attention to transformations between fixed and 
moving reference frames. In Section IIIII the form of the 
three lowest invariants of KdV equations is derived for 
different forms of the equations. In Section m we show 
that the energy calculated from the definition H = T +V 
has no invariant form. Section El describes the varia¬ 
tional approach in a potential formulation which gives a 
proper KdV equation but fails in obtaining second order 
KdV equations. In the next section the proper invari¬ 
ants are obtained from Luke’s Lagrangian density. Sec¬ 
tion |VlI] summarizes conclusions on the energy for KdV 
equation. In section IVIIII we apply the same formalism 
to calculate energy for waves governed by the extended 
KdV equation (second order). We found that energy is 








2 


not conserved neither in fixed coordinate system nor in 
the moving frame. 

II. THE EXTENDED KDV EQUATION 

The geometry of shallow water waves is presented in 

Fig- IH 

In [1,, 21 we derived an equation, second order in small 
parameters, in the fixed reference system and with scaled 
nondimensional variables containing terms for bottom 
fluctuations. They will not be considered here. 


bottom 
undisturbed surface 



FIG. 1. Schematic view of the geometry. 


The explicit form of the scaling leading to equations © 

© is given by (l29l) . 

Problems with mass, momentum and energy conser¬ 
vation in the KdV equation were discussed in [2Cj re¬ 
cently. In this paper the authors considered the KdV 
equations in the original dimensional variables. Then 
the KdV equatios are 

3 c ch 2 „ 

Vt +CVx + oTOTx + -z-Vxxx = 0, (5) 

2 h o 

in a fixed frame of reference and 

3 c ch 2 . . 

Vt + + ~x~Vxxx = 0, (6) 

2 h 6 

in a moving frame. In both, c = \fgh 1 and © is obtained 
from © by setting x' = x — ct and dropping the prime 
sign. 

In our present paper we discuss the energy formulas ob¬ 
tained both in fixed and moving frames of reference for 
KdV ©,©,©,©■ There seem to be some contradic¬ 
tions in the literature because the form of some invariants 
and the energy formulas differ in different sources because 
of using different reference frames and/or different scal¬ 
ings. In this paper we address this problems. 

The second goal is to present some invariants for a 
KdV type equation of the second order ©■ 


For a flat bottom that equation reduces to the second 
order KdV type equation, identical with [4j, Eq. (21)] for 
P = a, that is, 


3 1 2 

Vt + Vx + a -rjr] x + P-rj 3x + a [--ifv* 


(1) 


/ 23 


+a/3 —r] x r] 2 X + + P —— »? 5 cc = 0. 


V 24 


12 


19 


360 


Subscripts denote partial differentiation. Small parame¬ 
ters a, P are defined by ratios of the wave amplitude a, 
the average water depth h and mean wavelength l 



Equation © was earlier derived in Q and called "the 
extended KdV equation". 

Limitation to the first order in small parameters yields 
the KdV equation in a fixed coordinate system 


III. INVARIANTS OF KDV TYPE EQUATIONS 


What invariants can be attributed to equations Q 
© and © © ? 

It is well known, see, e.g. [7j, Ch. 5], that an equation 
of the form 


dT dX _ n 
~dt + ~dx ~ ’ 


( 7 ) 


where neither T (an analog to density) nor A' (an analog 
to flux) contain partial derivatives with respect to t, 
corresponds to some conservation law. It can be applied, 
in particular, to KdV equations (where there exist an 
infinite number of such conservation laws) and to the 
equations of KdV type like ©. Functions T and X 
may depend on x, t, 77 , r] x , r] 2x , ■ ■ ■ , h, h x ,... , but not rj t . 
If both functions T and X x are integrable on (— 00 , 00 ) 
and lim X = const (soliton solutions), then integration 

x —>-ioo 

of equation © yields 


Vt+Vx + Ot ~Wlx + P = 0. 
Transformation to a moving frame in the form 
x = (x — t), t = t, fj = rj, 


d 


P> 5t 

since 

( 3 ) 



X x dx = X(oo, t) 


T dx = const. 


X{— 00 , t) = 0. 


( 8 ) 

( 9 ) 


allows us to remove the term r] x in the KdV equation in 
a frame moving withthe velocity of sound \fgh 

3 1 

Vi + a -fjfjx + P Qfj 3x = 0. (4) 


The same conclusion applies for periodic solutions 
(cnoidal waves), when in the integrals © , © limits of in¬ 
tegration (— 00 , 00 ) are replaced by (a, &), where b— a = A 
is the space period of the cnoidal wave (the wave length). 
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A. Invariants of the KdV equation 

For the KdV equation © the two first invariants can 
be obtained easily. Writing © in the form 


one immediately obtains the conservation of mass (vol¬ 
ume) law 


I*- 1 -* = / ndx = const. 


Similarly, multiplication of © by 77 gives 


IG" 2 ) + s G” 2+ G" 3 - G' 3 ” 2+ G’ m ”)=°’ 

( 12 ) 

(13) 


resulting in the invariant of the form 

/ OO 

rj 2 dx = const. 

-OO 


In the literature of the subject, see, e.g. 2fi| |. I is 
attributed to momentum conservation. 

Invariants /l 1 ),/! 2 ) have the same form for all KdV 
equations ©, ©, (IA2I) . ©, ©. 

Denote the left hand side of © by KDV(cc, t) and take 

3 t 7 2 x KDV(*,i) - l~r h x ^-KDV(o;,t). (14) 

3 a ox 

The result, after simplifications is 


9 / 3 1 /? 2 \ 9/9 4 1 2 

S ’ ~ la' 1 *) + teXS* 1 + 


(15) 


1 9 2 2 1 /3 2 1/3 


o ~z 1 ij 1 ^ 2 ^ 2 \ n 

-PVxV + V + TTT- V2x - 7T VxVSx ~ X~V x =0. 

18 a 9 a 3a ) 

Then the next invariant for KdV in the fixed reference 
frame © is 


r(3) 

1 fixed frame 


/ ( r?3 “ \a v y dx = const ‘ G 6 ) 


The same invariant is obtained for the KdV in the mov¬ 
ing frame ©. The same construction like OD but for 
equation © results in 

9 / 3 1/9 2 \ 9/9 4 1 2 

»G ‘5s’ , G + sG"' + v ,m (i7> 

2 ^ 3 ^ 1 /3 2 2 1 /3 2 G n 

—PVxV + V + T^ — V2x - n —VxV3x = 0. 

18 a 9 a / 

Then the next invariant for KdV equation in moving ref¬ 
erence frame 0 is 


y( 3 ) 

1 moving frame 


/ if ~\^ r, *) dX = con st - d8) 


The procedure similar to those described above leads 
to the same invariants for both equations © and © 
where KdV equations are written in dimensional vari¬ 
ables. To see this, it is enough to take 3r] 2 x kdv(x, t) — 


/ 3 2 1 \ 

6 

in + 4 a7 ? + g/3 Vxx) = o. 

(10) u 
9 


dt 


h 3 


d_ 

dx 


9c 

8 h 7 


'32 


(19) 


~ch 2 ml + ^ ch 2 r) 2 i) xx + ch 5 rf xx - ^ ch 5 r] x r) xxx ) = 0, 


( 11 ) whereas for equation © the flux term is different 


9 / 3 h 3 2^\ , 9 /9c 4 2 2 

my> - chm * 

+ ljch 2 m 2 i lxx + -^ch 5 ri 2 x - ^ch 5 rj x r] x 


( 20 ) 

= 0. 


But in both cases the same I^ invariant is obtained as 
/ ( r/3 ~ y 9 x) dx = const. ( 21 ) 


r(3) 

J dimensional 


Conclusion Invariants I^ have the same form for fixed 
and moving frames of reference when the transformation 
from fixed to moving frame scales x and t in the same 
way (e.g. x’ = x — t and t' = t). When scaling factors 
are different, like in m, then the form of in the 
moving frame differs from the form in the fixed frame, 
see Appendix [A] 

For those solutions of KdV which preserve their shapes 
during the motion, that is, for cnoidal solutions and single 
soliton solutions, integrals of any power of r](x, t) and any 
power of arbitrary derivative of the solution with respect 
to x are invariants. That is, 


/ OO 

(■ Vnx) a dx = const, 

-OO 


( 22 ) 


where n = 0 , 1 , 2 ,..., and a £ R is an arbitrary real 
number. Then an arbitrary linear combination of 
is an invariant, as well. 


B. Invariants of the second order equations 

Can we construct invariants for KdV type equations of 
the second order? Let us try to take T = rj for equation 
0 . Then we find that all terms, except i] t , can be 
written as X x , as 

3 a l 2 ( 3 2 
Vx + a-T)T) X + p-r) 3x + a I --rj t). 


3 

4 C 


/ 23 5 

\ o2 19 

1 —nxV2x + j^nmx 

) 360 ??5 G 

n 2 + jjfimx - \°?n 3 

C13 , 5 \ 

19 « 2 

360 G mx- 

(m"- + l2 m2 -) + 


dx 

( 23 ) 
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As (1251) depends on p and space derivatives and also 
since all those functions vanish when x —> ±oo, the 
conservation law for mass (volume) 


rj(x, t) dx = const. 


(24) 


holds for the second order equation. 

(Conservation law (El) holds for the equation with an 
uneven bottom, as well.) 

Until now we did not find any other invariants for the 
secod order equations. Moreover, we can show that the 
integral I (2 ) (fl3l) is no longer an invariant of the second 
order KdV equation ©• 

Upon multiplication of equation m by p one obtains 


0 = 


d 

'V 

, l_\ 

dt 

^ J 

dx 


V + ^atif + - F /d(--vl + VV2. 


32 


360 


1 


V + —rP \7iVxx-VxV3x+VV 4c 


(25) 


+-^ a P rfm* 


l 


-a/3 irn x r] 2 x ■ 


The last term in (El) can not be expressed as 
J^-X(p, 7j x , ...). Therefore p 2 dx is not a conserved 

quantity. 


IV. ENERGY 


The invariant I^ is usually referred to as the energy 
invariant. Is this really the case? 


A. Energy in a fixed frame as calculated from the 
definition 


The hydrodynamic equations for an incompressible, in- 
viscid fluid, in irrotational motion and under gravity in 
a fixed frame of reference, lead to a KdV equation of the 
form 


3 1 

Vi + Vx + a -rjfji + /3 -p 3x = 0. (26) 

We will find the function 

fx = V~ ^ap 2 + ^PVxx, (27) 

obtained as a byproduct in derivation of KdV, useful in 
what follows. For more details see Appendix [B] or [mI . 
Chapter 5]. Tildas denote scaled dimensionless quanti¬ 
ties. 

Now construct the total energy of the fluid from the 
definition. 


The total energy is the sum of potential and kinetic en¬ 
ergy. In our two-dimensional system the energy in origi¬ 
nal (dimensional coordinates) is 


E=T+V= 




dx 


(28) 


+ 




dx. 


Equations (El) and (H71) are obtained after scaling 0, 
a a. We now have dimesionless variables, according to 



x 





t = 


and 


V = pgh 2 l 



r* 1+af/ 


pydydx, 


i/Vgh ’ 

(29) 


(30) 


T = i pghH 



f* \-\~Oil) 



Note, that the factor in front of the integrals has the 
dimension of energy. 

In the following, we omit signs ~, having in mind that 
we are working in dimensionless variables. Integration in 
(l30l) with respect to y yields 


V = ^ gh 2 lp 
= \gh 2 ip 


( a 2 rf + 2 ay + l) dx 



(a 2 p 2 + 2ap) dx + 



(32) 


After renormalization (substruction of constant term 
dx) one obtains 


1 „ r°° 

V = -gh 2 lp / (a 2 ? 7 2 + 2 ap) dx. (33) 

" J — OO 


In kinetic energy we use the velocity potential ex¬ 
pressed in the lowest (first) order 

4>x = fx — ^Py 2 fXXX and (j) y = -pyf xx , (34) 

where f x was defined in (1771) . Now the bracket in the 
integral (1711) is 


^a 2 (/) x 2 + —— a 2 (/ 2 + Py 2 {—f x fxxx + fix)) ■ 

(35) 

Inegration with respect to the vertical corrdinate y gives, 
up to the same order, 


i r+°° 

T = 7^P9h 2 lJ a 2 [f 2 ( 1 + ap) 


+P{~fxfxxx + + ar]) 3 dx (36) 


/• + 00 


= 2 P9 h ~ l / a 


fx + a fx 7 l + (fxx — fxfxxx) 


dx. 
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In order to express energy through the elevation funcion 
?7 we use m- We then substitute f x = rj in terms of 
the third order and f 2 =g 2 ^ \otrj 3 + ^/3r]r] xx in terms 
of the second order 


T = —pgh 2 l 


= | pgh 2 l a 2 


r+°° I" / 1 „ 2 

2 ' n - -ag' 3 + ~Prm xx 

+a?7 3 + -ft {rf x - T]r] xx ) 

f+ " (p + ‘ V) * 

r+°° i 

+ ~P {vl + VVxx) dx 

J — OO d 


dx 


(37) 


The last term vanishes as 


/ +oo /*+oo /*+oo 

{vl+VVxx)dx= r%dx + rm x \±™- rgdx = 0. 

-oo J— oo J—oo 

(38) 

Therefore the total energy in the fixed frame is given by 

E to t=T + V = pg/i 2 / J (arj + (ag) 2 + ar?) 3 ^ dx 

=pg/i 2 Z^a/( 1 )+a 2 /^ + ^a 2 l4 3 ) + -^a 2 /3 J g^dxj 


The energy (39) in a fixed frame of reference has non in¬ 
variant form. The last term in (39) results in small de¬ 
viations from energy conservation only when l] x changes 
in time in soliton’s reference frame, what occurs only dur¬ 
ing soliton collision. This deviations are discussed and 
illustrated in Section VI E. 

The result (1391) gives the energy in powers of 77 only. 
The same structure of powers in 77 was obtained by the 
authors of a, who work in dimensional KdV equations 
© and ©. On page 122 they present a non-dimensional 
energy density E in a frame moving with the velocity U. 
Then, if U = 0 is set, the energy density in a fixed frame 
is proportional to arj + a 2 g 2 as the formula is obtained 
up to second order in a. However, the third order term is 
ja 3 g 3 , so the formula up to the third order in a becomes 


E 


ag + a 2 g 2 + -a 3 rj 3 . 


(39) 


This energy density contains the same terms like (l3iH) and 
does not contain the term rj 2 , as well. 

Energy calculated from the definition does not contain 
a proper invariant of motion. 


by Ali and Kalisch [ 8 ,Sect. 3] working in such system one 
has to replace 4> x by the horizontal velocity in a moving 

frame, that is by = afj-\afi 2 +P fj xx -±. 

Then repeating the same steps as in the previous subsec¬ 
tion yields the energy expressed by invariants 


E tot =pgh 2 l 


= pgh 2 i 



^arj+^ag) 2 + ^a 3 \ 

/(i) + l a 2 /( 2 ) + 1 

4 2 


3 1/3-2 

V ~o~V. 
3 a 



dx 

(40) 


The crucial term —^oPflf] 2 in (HU appears due to inte¬ 
gration over vertical variable of the term —fj xx supplied 

by (h ~ if- 


V. VARIATIONAL APPROACH 

A. Lagrangian approach, potential formulation 

Some attempts at the variational approach to shallow 
water problems are summarized in G.B. Whitham’s book 
(3|ED, Sect 16.14], 

For KdV as it stands, we can not write a variational 
principle directly. It is necessary to introduce a velocity 
potential. The simplest choice is to take 77 = p x . Then 
equation (© in the fixed frame takes the form 

3 1 

^Pxt + ^Pxx ^Q^Px^Pxx + ~P^xxxx — 0 • ( 41 ) 

2 0 

The appropriate Lagrangian density is 

/* 1 1 2 cf 3 P 2 , . . 

Lfixed frame := ~~^PtPx — x Px ~ ~^P X + Pxx • (42) 

Indeed, the Euler-Lagrange equation obtained from La¬ 
grangian (H21) is just (ED. 

For our moving reference frame the substitution 77 = 
p x into ED gives 

3 1 

'pxt V ~^0-PxPxx T ~Pp X xxx = 0 • (43) 

So, the appropriate Lagrangian density is 

c 1 Cf o P 2 / . ., 

•^moving frame := ~~^PtPx ~ “j 'Px T J^Pxx ■ (44) 

Again, the Euler-Lagrange equation obtained from La¬ 
grangian Hill is just El- 


B. Energy in a moving frame 


B. Hamiltonians for KdV equations in the 
potential formulation 


Now consider the total energy according to (28) calcu¬ 
lated in a frame moving with the velocity of sound c = 
\fgh. Using the same scaling (29) to dimensionless vari¬ 
ables we note that in these variables c = 1. As pointed 


The Hamiltonian for the KdV equation in a fixed frame 
© can be obtained in the following way. Defining gen- 

dC _ 

eralized momentum n = ——, where C is given by (14211 . 

opt 
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one obtains 


H = 


/ oo pC 

[ 7 up — C]dx = / 

-OO J — ( 


' —OO 
poo 


1 2 

. 1 

( 3 ? 2 \] 

b 77 

+i a 



2^ + 4 % 


dx. 


12 


2 

^Pxx 


dx 

(45) 


The energy is expressed by invariants fa 2 \ j( 3 ) so it is a 
constant of motion. 

The same procedure for KdV in a moving frame 0 
gives 



The Hamiltonian (HtH) only consists I^. 

The constant of motion in a moving frame is 


E = - 1 ^ = const. 
4 


(47) 


The potential formulation of the Lagrangian, described 
above, is succesful for deriving KdV equations both for 
fixed and moving reference frames. It fails, however, for 
the second order KdV equation dTJ). We proved that there 
exists a nonlinear expression of C(ipt, <f x , Vxxi ■ • ■)> such 
that the resulting Euler-Lagrange equation differs very 
little from equation 0- The difference lies only in the 
value of one of the coefficients in the second order term 
a/3 (H VxVix + -p^vr/sx)- Particular values of coefficients 
in this term also cause the lack of the I^ invariant for 
second order KdV equation, (see (HOO¬ 


VE LUKE’S LAGRANGIAN AND KDV 
ENERGY 

The full set of Euler equations for the shallow water 
problem, as well as KdV equations ©, (IA2I) . and sec¬ 
ond order KdV equation 0 can be derived from Luke’s 
Lagrangian [2Sj, see, e.g. js]. Luke points out, how¬ 
ever, that his Lagrangian is not equal to the difference of 
kinetic and potential energy. Euler-Lagrange equations 
obtained from L = T —V do not have the proper form at 
the boundary. Instead, Luke’s Lagrangian is the sum of 
kinetic and potential energy suplemented by the 4>t term 
which is necessary in dynamical boundary condition. 


A. Derivation of KdV energy from the original 
Euler equations according to f24j] 

In Chapter 5.2 of the Infeld and Rowlands book the 
authors present a derivation of the KdV equation from 
the Euler (hydrodynamic) equations using a single small 
parameter e. Moreover, they show that the same method 
allows us to derive the Kadomtsev-Petviashvili (KP) 


equation [29] for water waves [3(3, [30, HH, [34] and also 
nonlinear equations for ion acoustic waves in a plasma 
sin. The authors first derive equations of motion, 
then construct a Lagrangian and look for constants of 
motion. For the purpose of this paper and for comparison 
to results obtained in the next subsections it is convenient 
to present their results starting from Luke’s Lagrangian 
density. That density can be written as (here g = 1) 


i-i+i) ■ 


L = 


fo + + <t> 2 .) + ^ 


dz. 


(48) 


In Chapter 5.2.1 of 0 the authors introduce scaled 
variables in a movimg frame {e plays a role of small pa¬ 
rameter and if e = a = /?, then KdV equation is ob¬ 
tained). Then (for details, see Appendix iBl or 0, Chap¬ 
ter 5.2]) 


fa = 4>x = eft - e 2 y /«««, 

fa = —£f$ + e 2 ^f T + — - e 3 — f^ T . (49) 


Substitution of the above formulas into the expression [ ] 
under the integral in (1481) gives 

[ ] = Z ~ (^fr + 2 ft + ( 50 ) 

+ £3 y [ - /f?T + (fh ~ /?/??«)] + ^(e 4 ). 

Remark The full Lagrangian is obtained by integration 
of the Lagrangian density (l48ll with respect to x. In¬ 
tegration limits are (— 00 , 00 ) for a soliton solutions, or 
[a, b\ , where b — a = A-wave length (space period) for 
cnoidal solutions. Integration by parts and properties of 
the solutions at the limits, see allow us to use the 
equivalence J^(/ £ 2 £ - /c/«f)d£ = fZo 2 /«^- 
Therefore 

[] = z - £ fz + e 2 ^/ T + 2-/f + ( 51 ) 

+ e3 y [ _ /?u- + 2 /f{] + 0 ( e 4 ). 

Integration over y gives (note that 1 + rj =>• 1 + erf) 


L= -(1 + erf) 2 + (1 + erf) 


-eh + e 2 [f T + -JI 


+ ^(1 +£7?) 3 


^£ 2 /s« - ^hir + e 3 /| 5 


2 2 
Write (15*21) up to third order in e 

L = L (0) + ei (1) + £ 2 L {2) + £ 3 L (3) + 0(£ 4 ). 


(52) 


It is easy to show, that 

£ (0) = y i (1) =^-/u 

i (2) = fr + -rj 2 - T/ft- + -f 2 + (53) 

L (3) = nfr + ^nfl + 
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The Hamiltonian density reads as 
H = f r^f + ~ L 

dfr dfe T 


(54) 


\ + e (v~h) + e 2 (- rifu + 77/f + Ifstt 


+£3 (^/l + ^/«« + ^4) 


Dropping the constant term one obtains the total energy 
as 


£ = 


£ (V- fr) + e 2 ( \il 2 - 7 7 / € + 


-3 I j-2 1 


+£ tt’7/ € + 




(55) 


Now, we need to express and its derivatives by 77 
and its derivatives. We use m replacing a and /? by e, 
that is, 


, 1 2,1 

fi=V-^V +3^- 


(56) 


Then the total energy in a moving frame is expressed 
in terms of the second and the third invariants 


£ = - 


1 r°° 

£ 2 - y 2 dx + s 6 


,1 


V - o 7 ?? dx 


(57) 


Note that the term occuring in the third order 

invariant originates from three terms appearing in cj) 2 , 
ft 2 and (f>t (see terms /« and /^ in (HiUl l. 


B. Luke’s Lagrangian 


The original Lagrangian density in Luke’s paper (23} is 


fh(x) 

L= / p 

Jo 


fit + \(fil + fil) + gy 


dy. 


(58) 


After scaling as in els a 


h,~x~ri~y~ 

7 — 7 =F</>, x = Ti rj = -, y =-, t = 

tay/gh l a ri 


we obtain 

fit = gha fit , fil = gha 2 $1, 


i/Vgh ’ 

(59) 


b y =g h —fiy. 


( 6 °) 

The Lagrangian density in scaled variables becomes 
(dy = hdy) 


L = pgha 


1 


r»l+a ?7 


1 


dy 


+ ^Pgh 2 (l+ag) 2 . 


(61) 


So, in dimensionless quatities 

dy+ \ay 2 , 

(62) 

where the constant term and the term proportional to 77 
in the expansion of (1 + < 377) 2 are omitted. The form (l62l) 
is identical with Eq. (2.9) in Marchant & Smyth Q. 

The full Lagrangian is obtained by integration over x. 
In dimensionless variables (dx = l dx) it gives 



C = Eq 



1 

2 



dy+-ay 2 


dx. 


(63) 

The factor in front of the integral, E 0 = pghal = pghfil a, 
has the dimension of energy. 

Next, the signs ( ~ ) will be omitted, but we have to 
remember that we are working in scaled dimensionless 
variables in a fixed reference frame. 


C. Energy in the fixed reference frame 


Express the Lagrangian density by 77 and / = . 

Now, up to first order in small parameters 


= / - \fiy 2 fxx, 


■H ft 2 fiy“ fxxti 


= fx - 77 Py 2 fxxx, fiy = -fiyfxx ■ (64) 


1 

2 PV Jxxxi tyy 

Then the expression under the integral in (lC>21) becomes 


[] = /*— \fiy 2 fxxt + ^aff + ^a/3y 2 (-f x f xxx + f xx ) ■ 

(65) 

From properties of solutions at the limits 
{-fxfxxx + fL) => 2 fix- Integration of {65]) over y 
yields 


pgha 


1 


1 


1 


= [ ft + -zO-fx (1 + ar ]) ~ - fifxxt oC 1 + ar i) 


( 66 ) 


+ a / 3 ff x 3(1 + ar l ) 3 + 2 m l 2 ' 

The dimensionless Hamiltonian density is 
(ft§E + fxxt^- - L) 


H 1 1 1 

— = ^ a fx( l + ar l) + a Pfxx 3(1 + arj) 3 + -ay 2 

(67) 

Again, we need to express the Hamiltonian by 77 and its 
derivatives, only. Inserting 


fx = h — ^arj 2 + ^/ 3 y xx 


( 68 ) 


into (1671) and leaving terms up to third order one obtains 

H 


pgh 2 l 


ay 2 + + ^a(3( y 2 x + yy xx ) 


■ (69) 






























The energy is 



since the integral of the a/3 term vanishes. Here, in the 
same way as in calculations of energy directly from the 
definition (l39l) . the energy is expressed by integrals of ry 2 
and rf. The term proportional to ary is not present in 
GOD, because it was dropped earlier |3j. 


D. Energy in a moving frame 

Transforming into the moving frame 
x = x-t, i = at, d x =d s , d t = -d. J: + ady (71) 


— f f xx ’ /c fxXXl 4 *y — Pyfxx, 

(72) 

= -fx + \fiy 2 fxxx + a(/t - \py 2 fxxi)- (73) 


8 

Finally the energy is given by 


E 

pgh 2 l 


1 1 

= a 2 - ifdx + a 3 - 
4 7_oo 2 




1/9 2 

3 a 


dx 


, (79) 

since integrals from terms with /3, /3 2 vanish at inte¬ 
gration limits, and — j^ggxx =>• y^ry 2 . The invariant 
term proportional to ary is not present in (17911 . because 
it was dropped in (1621) . If we include that term, the 
total energy is a linear combination of all three lowest 
invariants, j( 1 ),/( 3 ),/( 3 ). 


Comment An almost identical formula for the energy 
in a moving frame, for KdV expressed in dimensional 
variables ©> was obtained in [20], That energy is ex¬ 
pressed by all three lowest order invariants 

i r°° i r 2 r°° „ 

£ = --c 2 pdx+-— / ry 2 dx (80) 

^ J — OO 4 /l J — oo 



as well. Translation of 1801) to nondimensional variables 
yields 

£g = ggh 2 l ^a 2 I < - 2 '> + ia 3 /( 3 )^ . 


Up to second order 


\ (a<j% + = \ [a/J + u/ 3 y 2 {-fxfxxx + fix)] 


= + a/3ry fxx- 


(74) 


H 


pgh 2 l 


~fx + 2 a/IJ (! + ar l) 

+ i(l + ary) 3 Q/3/ 2$s + a/3+ ^ary 2 


(77) 


Expressing f x by (1551) one obtains 
H 


pgh 2 l 


- 2 + T^fiVxx - ^a 2 ry 3 


+ a/3 ( -\yl - j^yy^j - 


(78) 


E. How strongly is energy conservation violated? 


Therefore the expression under the integral in (l62l) is 

[ ] = -fx+^Py 2 fxxx+a{fi-^Py 2 fxxi)+\ufl+oiPy 2 flx- 

(75) 

Integration yields 

= (~fx + aft + ^a/ 2 ^ (1 + ary) (76) 

+ i(l + ary) 3 0/3(/*** - f X xt) + a/3/ 2 0 -0ary 2 . 
Like in (1571) above, the Hamiltonian density is 



FIG. 2. Precision of energy conservation for 3-soliton solution. 
Energies are plotted as open circles (Ei) and open squares 
(E 2 ) for 40 time instants. 

The total energy in the fixed frame is given by equa¬ 
tion (1391) . Taking into account its non-dimensional part 
we may write 


Ei(t) = 


T+V 


ag + (ary) 2 + -(ary) 3 


ggh 2 l 

1 r°° 

= + a 2 1^ + — / (ary) 3 da; 

^ J —OO 


dx 


(81) 
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In order to see how much the changes of Ei violate energy 
conservation we will compare it to the same formula but 
expressed by invariants 

E 2 (t) = al w + a 2 J (2) + i« 3 / (3) . (82) 

The time dependence of E\ and E 2 is presented in 
Fig. [2] for a 3-soliton solution of KdV ©• Presented is 
time evolution in the interval t £ [—12,0]. The shape 
of the 3-soliton solution is presented only for three times 
t = —12,—6,0 in order to show shapes changing during 
the collision. 

For presentation the example of a 3-soliton solution 
with amplitudes equal 1,5, 1 and 0.5 was chosen. In Fig. [3] 
the positions of solutions at given times were artificially 
shifted to set them closer to each other. The plots in 
Figs. [5] and [3] for t > 0 are symmetric to those which are 
shown in the figures. 

For this example the relative discrepancy of the enregy 
Ei from the constant value, is very small 


8E = 


E^t =-12) - Ei(t = 0) 
Ei(t = -12) 


0.000258. 


(83) 


However, the E 2 energy is conserved with numerical pre¬ 
cision of thirteen decimal digits in this example. In a 
similar example with a 2-soliton solution (with apmli- 
tudes 1 and 0.5) the relative error (1551) was even smaller, 
with the value 5E ss 0.00014. This suggests that the de¬ 
gree of nonconservation of energy increases with n, where 
n is the number of solitons in the solution. 



FIG. 3. Shape evolutiom of 3-soliton solution during collision. 


VII. CONCLUSIONS FOR KDV EQUATION 

The main conclusions can be formulated as follows 

• The invariants of KdV in fixed and moving frames 
have the same form. (Of course when we have the 
same scaling factor for x and t in the transforma¬ 
tion between frames). 


• We confirmed some known facts. Firstly, that the 
usual form of the energy H = T + V is not always 
expressed by invariants only. The reason lies in the 
fact, as pointed out by Luke in [23], that the Euler- 
Lagrange equations obtained from the Lagrangian 
L = T— V do not supply the right boundary condi¬ 
tions. Secondly, the variational approach based on 
Luke’s Lagrangian density provides the right Euler 
equations at the boundary and allows for a deriva¬ 
tion higher order KdV equations. 

• In the frame moving with the velocity of sound 
all energy components are expressed by invariants. 
Energy is conserved. 

• Numerical calculations confirm that invariants 

J( 3 ) in the forms (HU), jT5]), (HD, (T5D 
are exact constants of motion for two- and three- 
soliton solutions, both for fixed and moving coordi¬ 
nate systems. In all performed tests the invariants 
were exact up to fourteen digits in double precision 
calculations. 

• For the extended KdV equation © we have only 
found one invariant of motion 1^ (1551) . 

• The total energy in the fixed coordinate system 
as calculated in (1551) is not exactly conserved but 
only altered during collisions, even then by minute 
quantities (an order of magnitude smaller than ex¬ 
pected). Details in figure caption of figure [5] 

A summary of these conclusions can be found in Ta- 


VIII. EXTENDED KDV EQUATION 

In this section we calculate energy formula corespond¬ 
ing to a wave motion governed by second order equations 
in scaled variables, that is the equation (U for the fixed 
coordinate system and the correponding equation for a 
moving coordinate system. As previosly we compare en¬ 
ergies calculated from the definition with those Luke’s 
Lagrangian. 


A. Energy in a fixed frame calculated from 
definition 

Now , instead of © we consider the second order KdV 
equation, that is (U called by Marchant & Smyth 0 
"extended KdV". 

In section HVl A. total energy of the wave governed by 
KdV equation, that is the equation © with terms only 
up to first order in small parameters was obtained in (1391) . 
In calculation according to eq. © the potential energy 
is expressed by the same formula (1321) as previously for 
KdV equation. In the expression for kinetic energy the 
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/ oo 

rf dx. * Formulas in this column are written in 

-OO 



Euler 

Luke’s 

Integrals 

Potential 

KdV 


equations 

Lagrangian 

T+V 

Lagrangian 

dimensional ' 

Fixed frame 

a/«+a 2 /( 2 > + la 3 ^(3) 

09]) 

a 2 lM + ho?yW 

(f70[l 

a7 (l) +Q 2 / (2) + l Q 3 r? (3) 

Hi 

1 7 (2)+1 q7 (3) 

S3 

a7 C 1 ) +Q 2 7 ( 2 ) + l a 3 J? (3) 

S3 

Moving frame 

Wi^+Wi^ 

l Q 2/(2) + l a 3j(3) 

_l a7 (l) + l a 2 7 (2) + l Q 3 7 (3) 

ia/ (3) 

_1 q7 (1) + 1 q7 (2)+1 q 2 7 (3) 


m 

GSD 

(SSI 

SSI 

SSI 


velocity potential has to be expanded to second order in 
small parameters 


= / - V 2 fxx + ^P 2 y 4 fx 


with derivatives 


<t>x = fx - \fiy 2 ixxx + iiP 2 y A f x 

9-V = Pyfxx + «/3 y^fxxxx- 


(84) 


(85) 


Integrating over y and retaining terms up to fourth order 
yields 


T = \pyh 2 i 


p + OO 


a 


fx + a7 lfx + (fxx ~ fxfxxx) 


+a/3(y fix ~ V fxfxxx) 


4“/3 ( 2Q f XXX ^^fxxfxxxx^~ g ^fxfxxxxx 


( 86 ) 

dx. 


Expression (1551) limited to first line gives kinetic energy 
for KdV equation, see (1551) . 

Now, we use the expression for f x (and its derivatives) 
up to second order, see e.g. B Eq. (2.7)], @, Eq. (17)] 


fx=y - + i/3 ri xx + 


(87) 


+«/3 ( jqV? 1 1 


mx 


l 


-J' Vx 


2 ' ) 10 
Insertion (1571) and its derivatives into (155| gives 


T = l -pgh 2 l 


P + OO 


y 2 + i«»? 3 + i/3 (y l + yyxx) 


3 2 4 0 /29 2 3 2 

a 7? + a/3 1-^ + -y y xx 

o2 ( 1 2 7 19 

3"/3 I 2 q^ xx 45 yxVxxx + ^^yyxx: 


( 88 ) 


dx. 


From properties of solutions at x —> ±oo terms with /3 
and /3 2 in square bracket vanish and the term with a/3 
can be written form. Finally one obtains 


T = i pgh 2 l j 


p + OO 


y 2 + —ay a — ^-a 2 y A — f—a(3yy. 


16 


24 


dx. 

( 89 ) 


Then total energy is the sum of (1551) and (1891) 


Etot = pgh l 


ay + (ay) 2 + -(ay) 3 (90) 




dx. 


The first three terms are identical as in KdV energy 
formula (1551) . the last two terms are new for extended 
KdV equation 0. 


B. Energy in a fixed frame calculated from Luke’s 
Lagrangian 

Calculate energy in the same way as in Section fVll C, 
but in one order higher. In scaled coordinates Lagrangian 
density is expressed by (1621) (here we keep infinite con¬ 
stant term) 


fl+ari 

L = pgh z l { I a 
1 


4>t 


l 




dy 


From 


+ 7/(1 + ay) 2 


we have 

(ft — ft ~Py fxxt H - 24^ 2/ fxxxxt • 


(91) 


(92) 


Inserting (l92l) and (1551) into (l9lT) . integrating over y and 
retaining terms up to third order one obtains (constant 
term -i id dropped) 

^m =a {^ + ft) + a (l '" 2 + r]ft + \%) ~ 

1 2 2 A 2 1 1 

~^ 2 a 3 ” ( g fxx ~~ 7 ^y fxxt ~ — fxfxxx 

T ! 2 q/^ fxxxxt T~a /3 (yfxx y fxxt y fxfxxx) 


+ a /3 ( fxxx 2Q fxxfx 


(93) 


1 r | 1 x x \ ^3 fxxxxxxt } 

~^yjxxxxt + YJq Jxjxxxxx J — P 5Q4Q f ’ 
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The the Hamiltonian density 

dL dL dL dL 

H — ft 757- + Jxxt 777--- L 


Then the expression under integral over y in (1311) becomes 
(in the following terms up to fourth order are kept) 

,,2 


ULj f ULj f ULj f ULj T 2 \ 

Wt + XXt dfxxt + f{4x)t 5 /( 4 .)* + f{6x)t df (6x)t ~ = 1 - 2 a/, + a 2 /, 2 + y 2 a 2 PfL (98) 


IS 


1 


+ a ?/3 g fxfxxx 


+a 3 /3 ( ~\nflx + \y fxfxxx 


P ( fxxx 3 ~ 2Q fxxfxxxx 


(94) 


-f- y aftfxxx y a (3fxfxxx fxxxxa 

4 2 /O 2 | 1 r2 1 /> n 1 « .» 

+y a p U Jxxx qJxxJXXXX 1 -^2 Jxjxxxxx 

After integration over y one obtains 
1 ,, / ,+0 ° 

T = ^Pdh 1 J [l + a {y - 2/,) + a 2 (~2yf x + f 2 ) 


Now, we use /,in the second order (1571) and its deriva¬ 
tives. Insertion these expressions into (1571) nd retention 
terms up to thired order yields 


H 


,2 2 1 „, 3„3 1 3 „ 4„4 


(95) 


~pghfi = —ay — ay - -a y + -a y 

2 o ( 1 2 1 

+a p I ■ --t] x - -yy x 

3 Q ( 29 2 3 2 

+ “ ^ V 48 “ 8 ?? % 

, 2a2 ^ 1 2 _ 7 19 

P 1 qqVxVxxx 360 y^xxxx 


The energy is obtained by integration of (1951) over x (us¬ 
ing integration by parts and properties of y and its deriva¬ 
tives at x —> Too). Then terms with a/3 and a/3 2 vanish. 
The final result is 


/• + 0O 


E = —pgh 2 1 


ay + {ay) 2 + ~{ay) 3 




(96) 


dx, 


the same as (l90l) but with the opposite sign. 


C. Energy in a moving frame from definition 

Let us follow arguments given by Ali and Kalisch [20l . 
Sec. 3] and used already in Section IIVI B. Working in 
a moving frame one has to replace (f> x by the horizontal 

velocity in a moving frame, that is, (j) x -• Then in a 

frame moving with the sound velocity we have 


4 >x = fx~ 5/3 y 2 fxxx + 77/ 3 2 y 4 f x 
(fry = / 3 yfxx T s /3 y 3 fxxxx- 


(97) 


1 


1 


T 0 a/3 fxxx fa yfx rri^P f xx 
3 60 

+a 2 /3 Q/ 2 , + V fxxx 2 fxfxxx 
+a /3 ( 77 / 3 . 3 . + 77 fxxx ~ Vfx fxxx) 
3 ~a (3 ( — fxxx ~ y^/xx/xxxx 


(99) 


-^2 V fxxxxx 


T rc ,fxfx 
60 


dx. 


Then insertion f x ([571) and its derivatives yields 


7*+00 


T = 7;P9 h 1 


— oo L 
2 


1 2 2 1 o 

-ay --ay- -a/3y x 




yrjx 


19 

180 


a$ 2 y x 


7 44 3 q f 7 2,^2 

—a r] +a 3 /3 ( — Wx 


+ a 2 /3 2 1 11 - 2 


30 


1 


+ ^a/3 y x 
36 


233 

360 


dx. 


y x y x 


y yxx 

119 
+ 360 


( 100 ) 


yy x 


where constant term is dropped. Using properties of so¬ 
lutions at x —> Too this expression can be simplified to 

^ + OOf 

9 9 a 3 3 1 

-a T] 


1 


T = -pgh lj 


—ay — -a 2 y 2 + -a 3 y 3 - —c 4 - 4 

1 2 1 4 ' 16 


+ ^u 2 /3y 2 x +-^a 3 [3yyl +^a 2 /3 2 y 2 xx dx. ( 101 ) 


Then total energy is 


/* + OO ‘ 


Etot = pgh 2 11 


1 


OO 7 A A 

-ay+-a y + -a 3 y 3 - —a y 


i 7 2/0 2 . 1 3 a 2,1 2 o 2 2 

+ PV X +^ a PWx P Vxx 


dx. ( 102 ) 


In special case a = j3 this formula simplifies to 


/* + <X> ‘ 


E tot = pgh 2 l 


' —OO 

+q4 ( 


1 1 2 2 

—ay H —ay 
2 4 


« 3 f ^ 3 + ^’ 7 2 


1 

24 


yy x 


1 

40 


Vxx 


dx. 


(103) 
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D. Energy in a moving frame from Luke’s 
Lagrangian 

Follow considerations in Section ED but with KdV 2 
equation ©• Transforming into the moving frame 
through (ITTI) we have now 

4> = f~ 7;Py 2 fxx + ^P 2 y 4 fxxxx, (104) 


— fx 2 fiy JXXX T 24^ y fXXXXX) 


by — Pljfxx g ft y fxxxxi 


4*t — fx fxxx 24^ y f'XXXXX 

+tt(/t - \py 2 fxxt T 2^"y f xxxx t')~ 


(105) 

(106) 
(107) 


Inserting (11041) (11071) into (l9lT) one obtains Lagrangian 
density in moving frame as (constant term 1 is dropped 
as previously) 

= a(,? - A) + “ 2 (? ?2 + * “ ” f ‘ + 


1 


6 


1 


+ ^a/3fxxx + a vfi+-gf s ~ T7^ a P f 


1 


120 


+ a 2 /? ( 6 fL 6 hxi + 2 dfxXX g hfxxx 

+ o; (3 Tjfxx — 2 Vfxxt d - 2 y fxxx — ^yfxfxxx 
OL P ^ 4Q fxxx T^fxxfxxxx d~ Y2Q fxxxxi 

24 V fXXXXX + Y2Q fX fXXXXCf^ • ( 108 ) 

Then Hamiltonian density 

dL dL dL 

H ft r\ r T fxxt o r “ T fxxxxi O r ” ^ (109) 

&Jt O J xxt OJxxxxt 

after insertion of (|108p into (|109p yields 

+ A) + “ 2 (-^ 2 + nh ~ \ 11 

fxxx 7 > ^ Vfx d" TTTTT^/d fxxxxx 

+ a ^ ^— — fxx — 2 yfxxx *T —fxfxxa^j (110) 

+ o; ^ ^ 2 yfxx 2 d fxxx T 2 Vfxfxxx 
d" a P 4 ^fxxx d~ fxxfxxxx 

4~ J xxxxx 1^0 f x f xxxxx 


In order to express (II 101) by 77 only we use fx in the form 
(El) and its derivatives. It gives 


H 1 2 2 1 a 3 3 3 

^m = -r v +e a to»-s aTI 


(in) 




Wx 


19 

360 


Vxxxx 


7 


H——a 4 77 4 — a 3 B ( — 7777 ? H- g 2 Vxx 

32 ' V 24 16 

2o2 /ll 2 233 119 

— ot p I —llxx 4" TjT^nxVxxx 4- 'Y^Wxxxx 


~ T2 a ^ rh 




Then energy is given by the integral 
/* + 00 


E = ggh z l 


+ \ a ^ xx ~ £ a V ( 112 ) 


+cV/3 


1 


( 5 2 
V48^ + A mxx 
7 


19 2 , 7 4 4 

+ 360 a/3 + 32 “ 77 


-« 3 /3 ( Zj-wl 4- Y^V 2 VxxJ - if2 a 0 3 Vx 


—a 


2«2 I H 2 
60 /xa 


233 


119 


Y20^ x ^ xxx 4 ~ j2Qyy xxxx 


dx. 


From properties of solution integrals of terms with 
ai/3, a/3 2 , cc/3 3 vanish and terms with a 2 /3, a 3 /3, a 2 /3 2 can 
be simplified. Finally, energy is given by the follwing 
expression 


E = ggh 2 l 


r +00 

' —CO 

7 


-^V-^aV + ^aV (113) 


dx. 


1 2 2 3 

73 3 

7 2! 

- a 77 — a 

4 

vs” + 

48 Vx ) 


~l8 a2/3r]x ~ ^ a3/3w ® _ 4^ a2 ^ 2?7 ®5 
In special case when /3 = a the result is 

<.+00 

E = ggh 2 l / 


4-a 4 ( —-? 7 4 - 7777 ?- dx. (114) 

\32 24 x 40 ,XX /J 

If the invariant term I ^ = J ag dx is dropped in (1901) 
or (19611 then the energy calculated in the moving frame 
(I113[l have the same value but with oposite sign. 


E. Numerical tests 

1. Fixed frame 

In order to check energy conservation for the extended 
KdV equation Q we performed several numerical tests. 
First, discuss energy conservation in a fixes frame. We 
calculated time evolution governed by the equation © 
of waves which initial shape was given by 1-, 2- and 3- 
soliton solutions of the KdV (first order) equations. For 
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FIG. 4. Example of time evolution of 3-soliton solution. 


presentation the following initial conditions were chosen. 
3-soliton solution have amplitudes 1.5, 1 and 0.25, 2- 
soliton solution have amplitudes 1 and 0.5 and 1-soliton 
solution the amplitude 1. The changes of energy pre¬ 
sented in Figs. [H] and [5] are qualitatively the same also for 
different amplitudes. An example of such time evolution 
for 3-soliton solution is presented in Fig. [4] 

Time range in Fig. []] contains initial shape of 3-soliton 
solution with almost separated solitons at t = 0, inter¬ 
mediate shapes and almost ideal overlap of solitons at 
t = 315. In order to do not obscure details the sub¬ 
sequent shapes are shifted verticaly with respect to the 
previous ones. Note additional slower waves after the 
main one which are generated by second order terms of 
the equation (HI), that were already discussed in @|. 



FIG. 5. Energy (non)conservation for the extended KdV 
equation in the fixed frame ©■ Symbols represent values 
of the total energy given by formulas or (1961) . Full square 

symbols represent the invariant /(% 

We see that the total energy for waves which move ac¬ 
cording to the extended KdV equation is not conserved. 
Although energy variations are generally small (in time 


range considered they do not extend 0.001%, 0.004% 
and 0.005% for 1-, 2-, 3-soliton waves, respectively) they 
increase with more complicated waves. For additional 
check of numerics the invariant I 1 ' 1 ' 1 = f_^ ai^(x,t)dx 
for the eaquation (©) was plotted as Mass. In spite of 
approximate integration the value of I H) was obtained 
constant up to 10 digits for all initial conditions. 


F. Moving frame 


Here we present variations of the energy calculated in 
a moving frame. The time evolution of the wave is given 
by the equation CD transformed with (ED, that is the 
equation 


3 1/3 

+ -7~V3x 
2 b a 


(115) 


23 


- -a 7 pfe + p —rjxmx + -rz’nmx + r/ 5 x = 0. 


24 


12 


19 p 2 


360 a 



FIG. 6. Energy (non)conservation for the extended KdV 
equation in the moving frame <fTT5l) . Symbols represent values 
of the total energy given by the formula (11021) . Full square 
symbols represent the invariant HV 

The time range of the evolution was chosen for a conve¬ 
nient comparison with the numerical results obtained in 
fixed reference frame, that is 2- and 3-soliton waves move 
from separate solitons to fully colliding time instant. The 
convention of symbols is the same as in Fig. [5] the energy 
is calculated according to the formula (11021) . In moving 
coordinate system energy variations are even greater than 
in the fixed reference frame, because in the time period 
considered it approaches values of 0.02%, 0.12% and 0.2% 
for 1-, 2- and 3-soliton waves, respectively. This increase 
of relative time variations of energy can not be atributed 
only to two times smaller leading term (i arp in (11021) 
with respect to (1901) . Again, in spite of approximate in¬ 
tegration the value of I^ was obtained constant up to 
10 digits for all initial conditions. 
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1. Conclusions for extended KdV equation 


We calculated energy of the fluid governed by the ex¬ 
tended KdV equation ID in two cases: 


1. In a fixed frame 


and I VIII B|) . 


2. In the frame moving with the sound velocity (sec¬ 
tions IVIII Cl and IVIIIDI) . 

In both cases we calculated energy using two methods: 
from definition and from Luke’s Lagrangian. Both meth¬ 
ods give consistent results. For fixed frame energies (1*1(1 
and (1961) are the same. For moving frame the energy 
calculated from the definition contains one term more 
then energy calculated from Luke’s Lagrangian, but this 
term (f arjdx) is the invariant /B). When this term is 
dropped both energies in moving coordinate system (11021) 
and (11131) are the same and energies in both coordinate 
systems differ only by sign. 

The general conclusion concerning energy conservation 
for shallow water wave problem can be formulated as fol¬ 
lows. Since there exists the Lagrangian of the system 
(Luke’s Lagrangian) then exact solutions of Euler equa¬ 
tions have to conserve energy. However, when approx¬ 
imate equations of different orders resulting from exact 
Euler equations are considered, energy conservation is 
not a priori determined. The KdV equations obtained in 
first order approximation has a miraculous property, an 
infinite number of invariants with energy among them. 
However, this astonishing property is lost in second or¬ 
der approximation to Euler equations and energy in this 
order may be conserved only approximately. 


Appendix A 

The simplest, mathematical form of the KdV equation 
is obtained from (ID by passing to the moving frame with 
additional scaling 

z=yfor-*), u = 7? ’ ( A1 ) 

which gives a standard, mathematical form of the KdV 
equation 

Ui A 6 u u s + —Uxxx = 0, or 
a 

uj + 6 u Ux + Uxxx = 0 for /3 = a. (A2) 

Equations (1A2|) , particularly with j3 = a are favored by 
mathematicians, see, e.g. [25J. This form of KdV is the 
most convenient for ISM (the Inverse Scattering Method, 
see, e.g. (2p(28j). 

For the moving reference frame, in which the KdV 
equation has a standard (mathematical) form (IA2I) . the 
invariant /is slightly different. To see this difference 
denote the lhs of (IA2I) by KDVm(a;, t ) and construct 

3?7 2 x KDVm(;r,f) — —r] x x ^-KDVm(a;, t) = 0. 
a ox 


Then after simplifications one obtains 


d 


1/3 


2 a 


d 


~-V x + yr- 


dx 


9 4 R P 2 
-r] -Q-wix 
2 a 


(A3) 


,o@ 2 1 (P 

+3 ~V dxx — r I -n. 
a 2 V a 




a 


T]x 1Jxxx 


= 0 , 


which implies the invariant I^ in the following form 


r(3) _ 

L moving frame 

r(3) 

1 moving frame 


= / ( rf — -—A?) dx = const or (A4) 

J -ooV 2a ) 


T] — -rj x ) dx = const for (3 = a. 


We see, however, that the difference between (1X41) and 
(1151) is caused by additional scaling. 

In the Lagrangian approach as described in Sect. El 
the substitution u = (p x into (IA2I) gives 


*Pxt A b (px^Pxx A <Pxxxx — 9 ■ (A5) 


Then the appropriate Lagrangian density for equation 
(1A2I) with (a = /3) is 

■^standard KdV := ~ Px A T^Pxx ■ (A6) 


Indeed, the Euler-Lagrange equation obtained from the 
Lagrangian (IA6I) is just (1A5I) . 

The Hamiltonian for KdV (1A2I) can be found e.g. in 


[22| . Defining generalized momentum n 
C is given by (IA6I) . one obtains 


dC 

dipt 


where 


H = 


r [ t ip - c\ dx = f 

J — oo J — 


dC r 

-oo \_Oipt 


Px .-^Pxx 


dx = 


1 


- ~ 2 r,l 


dx (A7) 
dx. 


This is the same invariant as 


j( 3) 

^moving frame 


in (IA4I) . 


Appendix B 


The set of Euler equations for irrotational motion of 
an incompresible and inviscid fluid can be written (ne¬ 
glecting surface tension) in dimensionless form: 


<1 

to 

-e- 

II 

o 


(Bl) 

<\>z= 0 

on z = 0 

(B2) 

Vt H - tfixVx 4*x = o 

on z = 1 + r] 

(B3) 

<!>t + \ (4>l + <t>V) + V = 0 

on z = 1 + rj. 

(B4) 

We look for solutions to the Laplace equation (1B1I) in 
the form 

OO 

<t>=Y'Z n f ln) (x 

n =0 

',y,t ) 

(B5) 
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yielding 

OO 

J2 [«(n - 1 )z n ~ 2 3 4 5 6 7 8 9 / (n) + z n V 2 f M 

n =0 

In two dimensions (x, z) we obtain 

f{n+2) _ 


= 0. 


-1 


Q2 j(n) 


(n + l)(n + 2)) dx 2 


(B6) 


(B7) 


The boundary condition at the bottom, <f> z = 0 at z = 0 
implies /B) = 0 and then all odd f( 2k + 1 ') = 0. Now 


= £(-!>’ 


J1 m f)2m 


d 2mf 


n —0 


(2m)! dx 2m ’ 


(B8) 


where / := f(°\ In the stretched coordinates d\ = ed\ 


so 


= ^1/ + £(-!)" 


n—0 


(2m)! 


K) 2m / 


(B9) 


Now both (1B1I) and (IB2I) are satisfied. We must also 
satisfy the boundary conditions on 2 = 1 + ?y. 

In the derivation of KdV and Kadomtsev-Petiashvili 
[29] from the Euler equations (IB 1 1) (IB4I) Infeld and Row¬ 
lands 0 applied scaling assuming the following relations 

vawelength : depth : amplitude as e ~ 1 ^ 2 : 1 : e. 


They then applied a transformation to a frame moving 
with velocity of sound. The coordinates scales as 


£ = £2 (x-t), 


T = £2 t 


(BIO) 


d t = — s 2 + e 2 d T , d x = e 2 d (Bll) 

For the wave amplitude and velocity potential the appro¬ 
priate scaling was 


ij = erf 1 ) + e 2 f 2 ^ + .... 


and 


(j) = £2 (ff) -+- £2 0 ( 2 ) 

Then the lowest order expression for cj> is 


(B12) 


(B13) 


<j>~£ 2 f- £5 (B14) 

Next, Infeld and Rowlands show that in order to simul¬ 
taneously satisfy (IB3I) and (IB4I) the next order contribu¬ 
tions to r] and (f> cancel. It is enough to keep 

1 T r] = 1 T and (j> = e (B15) 

and drop upper index B) i n what follows. 
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